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ABSTRACT 

We analyze results for the Boson Peak from the neutron time of flight spectroscopy data 
on Ge-As-Se, and Raman spectra data on m-TCP and OTP, using a recent mode coupling 
model that takes into account the coupling of density fluctuations with vibrational modes 
in presence of defects in the supercooled state. From the experimental results for different 
materials we observe that for more fragile systems characterized by increasing fragility 
parameter m, a slower relaxation of the defect-density correlation is needed to give rise 
to the observed peak in the spectra. 



PACS number(s) : 64.70P, 05.60, 64.60C 



The extra intensity observed in undercooled liquids in the neutron scattering [jl], Q as 

well as in Raman scattering p, Q at low frequencies, distinct from the quasi-elastic peak 

is usually referred to as the Boson peak in the literature. This characteristics feature of 

the supercooled liquid has been ascribed to the coupling between the relaxational and 

vibrational motions in the supercooled liquid in a recent work |^ with self consistent 

mode-coupling ^ model. The mode coupling theory has a better agreement with 

dynamics in fragile glassy systems. The classification ||^ used for glassy systems as 

strong and fragile depends on the ease with which structural degradation occur in those 

systems. A convenient measure for the fragility of the system is computation of the 

slope m of the relaxation time r curve against temperature T scaled with respect to the 

corresponding glass transition temperature Tg, 

dlog < r > 
d(T,/T) "^-^^ 

In Ref. IP we have described an extension of the simple mode coupling formalism to 
include the distinct vibrational modes that develop at low temperatures in the amor- 
phous state for understanding the extra intensity appearing for the structure factor. The 
model follows from the equations of Nonlinear Fluctuating Hydrodynamics ||^ extended 
to include the defect density in an amorphous solid like system which also sustains trans- 
verse sound modes. Here we have applied the model to extensive data comparison 
to understand the underlying relationship between the fragility of the system and the 
Boson Peak formation. We fit the model with respect to the data of Russina et. al. 
from neutron time-of-flight-spectroscopy for GeO.033 AsO.033 SeO.934 glass, as well as 
the data for m-tricresyl phosphate (m-TCP) by Sokolov et. al. from Raman scattering. 
The present analysis demonstrates that the criteria for the appearance of the peak is 
crucially related to the dynamics of defect densities in the disordered system. 

In studying the feedback effects on dynamics due to slowly decaying density fiuctua- 
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tions at supercooled states, the memory function H[(f){t)\ is obtained as a functional of 
the hydrodynamic correlation functions, in the following q independent form, 

H{t) = CiF[</)i(t), + C2^' W (2) 

where Ci and C2 are dimensionless constants determined in terms of the wave vector 
integrals due to the mode coupling vertex functions. 0^ and 0t are the correlation 
functions for the longitudinal and transverse sound modes. The function F{t) is expressed 

as m, 

F[m = e-'' + f{cr)Mt) (3) 

where 6 represents the time scale of very slowly decaying defect density and f{a) = 
(12 - 14a)/9(l - 2a) with a = {3\ - 2/i)/[2(3A + /i)] is the Poisson's ratio. For more 
details we refer to Ref. 0. Following the procedure described there we have used this 



model to fit the data of Russina et. al. |10, 11] for Ge-As-Se. The central quasi-elastic 



peak is fitted a Lorentizian of width 6, and is equivalent to taking the time scales of 



relaxation of the defects to be same [0, |T3| as that of the final time scale of relaxation of 
the density fluctuations. We show in Fig. la the dynamic structure factor as a function 
of frequency. The time scale for the decay of the defect density denoted by 6 plays a 
central role in appearance of the peak on the shoulder of the quasi-elastic peak. This 
intermediate peak disappear in the shoulder of the quasi elastic peak as 6 become large. 
This is shown by the fit of the data of Russina et. al. for a higher temperature in 
Figure lb. We analyze the data for temperatures, T = 252, 334, 359, 402, 440 and 502 
degree Kelvin, with suitable values of ci and C2 in the expression (0) for the memory 
function. The parameter 6 which gives the time scale of relaxation of defects is used 
here as the only adjustable parameter. It is considered in inverse units of the time Tq 
in terms of which the MCT equations are expressed in the dimension less form. Tq can 
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be obtained in terms of the microscopic frequencies of the hquid state. To study the 
behavior of 6 with the fragihty parameter we have considered data on another material 
m-TCP with a different m, used by Sokolov et. al. |Q and done a similar fitting to the 
Boson peak at temperatures T = 205, 235, 262 and 287 degree Kelvin. In Fig. 2 the 
fit to the Boson Peak data Q for m-TCP is shown. Variation of log(5) with the inverse 
temperature Tg/T is shown in Fig. 3, where Tg is the glass transition temperature of the 
corresponding material. The data points for m-TCP are shown by stars (*) and for OTP 
by filled circles (•). Variation of log{6) with temperature for Ge-As-Se alloy is shown 
in the inset by squares (□). Solid lines show the straight line fit to the data points. For 
lower temperatures the quantity 6 is small indicating that the defect densities are long 
lived and the Boson Peak appears to be more pronounced. The more fragile the system 
is, sharper is the fall of 6 which represents the time scale of defect correlation. In all three 
cases 6 shows an Arrhenius fall with temperature and a corresponding activation energy 
can be computed from the slope of the curve. In figure 4 we show the plot of activation 
energy A (in units of temperature T) for different systems with the corresponding fragility 
index m. With the systems of increasing fragility, more long lived defects are needed to 
give rise to the observed peak in the spectra. 

We have approximated here through S relaxation of the defect density by single ex- 
ponential mode. The full wave vector dependence has to be considered to account the 
coupling of the structural relaxation to the vibrational modes. The explicit temperature 
dependence of the peak is not captured in the present model. This can be computed 
through proper input for the static or thermodynamic properties that appear in the 
mode coupling integrals. Figure 3, demonstrates the key result of this paper that for 
more fragile systems 6, inverse of which relates to the defect density correlation, shows a 
sharper fall with temperature. Also we like to point out here that the temperature range 
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covered for the materials in this paper actually correspond to the part where the fragile 
glasses starts showing a sharp increase of viscosity on the Angell plot[^ - more fragile 
the liquid is, more dramatic is the increase giving a higher value for the fragility index 

m am. 
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Figure Captions 
Figure 1 



The neutron scattering data ( in arb. units) of Ref. [10, 11| normalized with respect to 
the Bose factor uj[n{uj) + 1] ( open circles) at (a) T = 252° K and (b) T = 440° K, Vs. 
the frequency in Mev. The solid line presents the result obtained from the present model 
for the normalized correlation function ip. 

Figure 2 

The Raman Spectra data ( in arb. units) for TCP normalized with respect to the Bose 
factor uj[n{uj) + 1] ( open circles) at T = 205° K, vs. the frequency in GHz. The solid 
line presents the result obtained from the present model for the normalized correlation 
function ip. 

Figure 3 

log(5) as a function of temperature [Tg/T) for OTP (•) and for m-TCP(*). Inset: For 
Ge-As-Se alloy (□). The solid lines represent the straight line fit to the data points. 

Figure 4 

Slope for the 5-1 /T curve, A ( in unit of "K ) Vs the fragility parameter m. 
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